ELb473fl7371US 

Adaptive Generalized Matched Filter Rake Receiver System And Method 

CROSS-REFERENCE TO RELATED APPLICATION 
This application claims priority from and is related to the following prior application: 
5 Adaptive Generalized Matched Filter Rake Receiver System And Method, United States 
Provisional AppUcation No. 60/257,737, filed December 22, 2000. This prior appUcation, 
including the entire written description and drawing figures, is hereby incorporated into the 
present application by reference. 

m BACKGROUND 
^ 1. FIELD OF THE INVENTION 

This invention relates generally to the field of spread spectrum rake receivers. More 
m particulariy, an Adaptive Generalized Matched Filter rake receiver system and method is 
J provided that is especially well suited for use in a mobile communication device. 

2. DESCRIPTION OF THE RELATED ART 

Mobile conmiunication devices operate in a multi-path propagation environment, i.e., 
there is typically more than one propagation path from the transmitter to the receiver. In 
addition, the velocity of the mobile device may vary from 0 km/h (standing still) to 500 km/h 
20 (traveling in a high speed train). Therefore, the multi-path propagation environment will 
typically range from direct line of sight to multi-clustered, multi-path propagation with no direct 
line of sight spread over several microseconds. Consequently, typical mobile communication 
devices employ a multi-fingered rake receiver that uses simple Maximal Ratio Combining and 
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standard pilot tracking processing in order to track the centroids of the multi-path clusters in a 
spread spectrum signal, such as a Code Division Multiple Access (CDMA) signal. 

A typical Maximal Ratio Combining (MRC) rake receiver includes a plurality of fingers, 
each of which correlates to a different delay of an input signal. The correlator outputs from each 
5 finger are then typically weighted by a vector of complex weighting coefficients, and combined 
to form a decision variable. In typical MRC rake receivers, the values of the coefficients in the 
weighting vector are chosen without regard to the statistical correlation properties of the noise 
impairment in the received signal, for instance by setting each weighting coefficient as the 
O complex conjugate of the channel impulse response. As a result, typical MRC rake receivers 
to perform optimally when the noise corruption to the input signal is limited to Independent 
Additive Noise (IAN), such as Additive White Gaussian Noise (AWGN), which is independent 
^ of the signal transmitted to the mobile device from a base station. In typical mobile 
fU concmiunication systems, however, multiple spread spectrum signals are transmitted at a single 

T, 

^ bandwidth, resulting in dependant noise, such as Multi-User Interference (MUI). Because 
tS typical MRC rake receivers are optimized to compensate for IAN, they are often sub-optimal 
when dependent noise is present. 

The use of a Generalized Matched Filter (GMF) to compensate for dependant noise in a 
spread spectrum signal is known. For instance, a generic description of a GMF is found in Kay, 
"Fundamentals of statistical signal processing - detection theory," Prentice Hall, 1998. In 
20 addition, the use of a GMF in a CDMA receiver is disclosed in G. Bottomly et al, "A generalized 
Rake receiver for interference suppression," IEEE Journal on selected areas in conmiunications, 
Vol. 18, No.8, Aug. 2000. Known Generalized Matched Filters, however, require an excessive 
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amount of processing, and are therefore not typically implemented in mobile communication 
devices. 

SUMMARY 

5 An Adaptive Generalized Matched Filter (AGMF) rake receiver system includes a rake 

receiver and an AGMF weight determination module. The rake receiver is coupled to a spread 
spectrum input signal and applies a vector of weight signals to the spread spectrum input signal 
to compensate for dependant noise and to generate a decision variable. The AGMF weight 
^ determination module monitors the decision variable and generates the vector of weight signals, 
IQ wherein optimal values for the vector of weight signals are calculated by the AGMF weight 
Ijl determination module by varying the vector of weight signals until the signal to noise ratio of the 
5 decision variable reaches a peak value. 

*r BRIEF DESCRIPTION OF THE DRAWINGS 

?s=s? 

Li' 

Fig. 1 is a block diagram of an exemplary rake receiver for use in an Adaptive 
Generalized Matched Filter (AGMF) rake receiver system; 

Fig. 2 is a block diagram of an exemplary AGMF rake receiver system; 

Fig. 3 is a graph plotting exemplary multi-path components {hi, hi, hj} of the channel 
impulse response h as a function of time; 
20 Fig. 4 is a more detailed block diagram of the AGMF Weight Determination module 

shown in Fig. 2; 

Fig. 5 is a graph plotting the SNR of the decision variable output z from the rake receiver 
as a function of the scalar parameter ro; 
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Fig. 6 is a block diagram of an exemplary Dual Decision Statistic Pilot Rake Receiver; 

and 

Fig. 7 is a flow diagram illustrating an exemplary method for calculating the optimal 
weight signal vector wopt . 

5 

DETAILED DESCRIPnON 
Referring now to the drawing figures. Fig. 1 is a block diagram of an exemplary rake 
receiver 10 for use in an Adaptive Generalized Matched Filter (AGMF) rake receiver system. 
O The rake receiver 10 includes a plurality of correlator fingers 12, a plurality of weight multipliers 

iB 14, and an adder 16. Each correlator finger 12 includes a delay element 18, a multiplier 20 and 

H 

Tz an integrator 22. Also shown is a receiver chain impulse response block 24, a pulse shaping 

y'l 

^' filter 26 and a coded sequence Co(n) block 28. 

fij A multi-path, spread spectrum signal x(t), such as a CDMA signal, is received by the 

£ mobile communication device, and is filtered by the receiver chain impulse response block hrx(t) 
i§ 24 to generate a demodulated base-band input signal r(t) to the rake receiver 10. The receiver 
chain impulse response block hrx(t) 24 represents the combined filter responses in the receiver 
chain prior to the rake receiver 10, such as the responses from an RF filter, band limiting 
components, an IF filter, and DSP filtering blocks. The input signal r(t) to the rake receiver 10 is 
coupled to one input of the multiplier 20 in each correlator finger 12. 
20 Each correlator finger 12 also receives a despreading signal c(t), which is formed by 

convolving the coded sequence Co(n) 28 for the desired traffic channel (n is the chip index) with 
the impulse response hp(t) of the pulse shaping filter 26. The impulse response hp(t) may be a 
single impulse or a rectangular pulse, depending upon how the correlation function is 



4 



implemented. Each correlator finger 12 is represented with an index number 1,2,,..., J. 
Operationally, each correlator finger 12 is substantially the same as the first, which will be 
described next in greater detail. A delay element (di) 18 is then applied to the despreading signal 
c(t) within each correlation finger 12 in order to generate a shifted despreading signal c(t-di) that 
5 is aligned with one channel of the multi-path input signal r(t). The shifted despreading signal 
c(t-di) is coupled to a second input of the multiplier 20. The multiplier 20 performs a complex 
operation on the input signal r(t) and the shifted despreading signal c(t-di), forming the product 
^ c(t-di)*r(t), where denotes the complex conjugate operation. The output of the multiplier 20 
p is then coupled to the integrator 22 in order to correlate the signals over some period of time and 
W to generate a correlation output y(di). The other correlator fingers 12 operate in substantially the 
III same way as described above, except that delay di is substituted with delay d2,...,dj for each of 
f the other correlator fingers 12. 

If the propagation channel of the input signal r(t) were ideal, i.e. a single-path 
p environment with no noise, then the rake receiver 10 would only require one correlation finger 
15 12 and a single delay element di. In this ideal case, the delay element di would be calculated 
such that the shifted despread signal c(t-di) would align exactly with a pilot signal within r(t), 
satisfying the equation: 
hic(t-di) = r(t), 

where hi, in this ideal case, is a single complex constant. Then, assuming that the correlation 
20 epoch is appropriately chosen, the correlation output y(di) reasonably approximates hi. Thus, 
the correlation output y(di) is an estimate of the channel impulse response of the complete link 
from the transmitter to the receiver. 
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In a true multi-path environment, however, the channel impulse response from the 
transmitter to the receiver is represented by a series of impulses of amplitudes {hi, ha, hj}, 
designated hereinafter by the vector h. Thus, a series of delays {di, da, dj}, represented 
hereinafter by the delay vector d, should be calculated for the array of correlation fingers 12, 

5 resulting in an array of correlator outputs {y(di), y(d2), y(dj)}. When the delay vector d is 
applied to a traffic-carrying input signal r(t), the array of correlator outputs may be approximated 
as follows: 

y(di) = Shi 

O y(d2) = Sh2 

t6 

g y(dj) = Shj, 

where S is an unknown complex amplitude coefficient that reflects the data content of the input 
m signal r(t). 

In order to estimate the value of the coefficient S, the correlator outputs are weighted by a 
k5 vector w of complex weight signals {wi, W2, wj} in the weight multipliers 14. The outputs 
from the weight multipliers 14 are combined in the adder 16 to generate a decision variable z 
which is proportional to the complex amplitude coefficient S by a real constant of 
proportionality. A system and method for deriving optimal values for the weight signals {wi, 
W2, wj} is discussed in detail below with reference to Figs. 2-6. 
20 Fig. 2 is a block diagram of an exemplary AGMF rake receiver system 30. The system 

30 includes the rake receiver 10, a CDMA processing module 32, an AGMF weight 
determination module 34, and a decoder 36. The CDMA processing module 32 and the AGMF 

weight determination module 34 respectively calculate the delay d and weight signal w vectors 
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used by the rake receiver 10. The CDMA processing module 32 and the AGMF weight 
determination module 34 are preferably software modules executing on a processing unit, such 
as a microprocessor, a field progranmiable gate array (FPGA), a digital signal processor, or a 
software interpreter module. It should be understood, however, that the exemplary AGMF rake 
5 receiver system 30 is not hmited to an embodiment having independent software modules for the 
CDMA processing module 32 and the AGMF weight determination module 34. Rather, the 
functions of the CDMA processing module 32 and the AGMF weight determination module 34 
may be performed by a plurality of separate software modules, by the same software module, or 
O by some other processing means. 

pb The CDMA processing module 32 receives the demodulated base-band input signal r(t) 

Jji as an input, and calculates the channel impulse response h and the delay vector d . The CDMA 
« processing module 32 tracks the CDMA forward-link pilot channel and measures the impulse 
ru response fi of the propagation channel. Contained in this impulse response h are the resolvable 
J multi-path clusters or components {hi, h2, hj} that are tracked in order to calculate the delay 
15 vector d , which is applied to the traffic-carrying input signal r(t) in the rake receiver 10. It 
should be understood, however, that alternative processing modules may be utilized in place of 
the CDMA processing module 32 that are configured for standards other than the CDMA 
standard. 

Fig. 3 is a graph 38 plotting exemplary multi-path components {hi, h2, hj} of the 
20 channel impulse response h as a function of time. This graph 38 illustrates that by tracking the 
multi-path components {hi, h2, hj} of the channel impulse response h, the CDMA processor 
module 32 may calculate the delays (di, d2, dj} between ±e multi-path clusters of the input 
signal r(t). 
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Referring again to Fig. 2, the AGMF weight determination module 34 receives the 
channel impulse response fi and the delay vector d from the CDMA processing module 32 and 
a feed-back decision variable z from the rake receiver 10, and generates the vector of weight 
signals w . A detailed description of the AGMF weight determination module 34 is provided 
below with reference to Figs. 4 and 5. 

An embodiment of the rake receiver 10 is described above with reference to Fig. 1. It 
should be understood, however, that the rake receiver 10 described above with reference to Fig. 1 
is just one possible embodiment, and may be replaced with the Dual Decision Statistic Pilot Rake 
Receiver 70 described below with reference to Fig 6, or with other rake receiver designs. 

The decision variable output z from the rake receiver 10 is also coupled as an input to the 
decoder 36, which converts the decision variable z into a binary receiver output Bout- The 
decoder 36 is preferably chosen based on the type of modulation scheme expected in the input 
signal r(t). For instance, the CDMA, DS-CDMA and UTMS standards typically employ 
quadrature amplitude modulation (QAM) schemes, while other standards, such as the GSM and 
GPRS standards, typically employ GMSK modulation. 

Fig. 4 is a more detailed block diagram of the AGMF Weight Determination module 34 
shown in Fig. 2. The AGMF Weight Determination module 34 preferably includes various sub- 
modules, including a total noise covariance matrix R„ module 40, a dependent noise covariance 

matrix Rdep module 42, an independent noise covariance matrix Rjp^ module 44, an optimizer 

module 46, a signal-to-noise ratio (SNR) module 48 and a weight determination module 50. 
These sub-modules may be independent software modules or sub-routines within the CDMA 
processing module 32, or may be realized by some altemative progranmiing structure or 
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processing device. The various sub-modules 40-50 are used by the AGMF Weight 
Determination module 34 to calculate the optimal weight signals Wopt . 

With reference to Fig. 1, the vector Y of rake correlator outputs {y(di), y(d2), y(dj)} 
may be modeled by the linear relation: 
5 Y=h+U 

where U is a noise vector. In order to achieve an optimal SNR for the decision variable z, the 
weight signals w should be calculated according to the following equation: 
Wopt = R„~^h 

p where R„ is the total noise covariance matrix for the noise vector U . 

f$ The noise vector tJ includes two relative components for the purposes of the AGMF 

Weight Determination module 34: an independent noise component, Uind, and a dependent 
N= noise or multi-user interference (MUI) component, Udep . Thus, the noise vector tJ may be 
H expressed as: 
^ U = Uind + Udep . 

15 The covariance matrix of tJ can be expressed by the superposition of the covariance 

matrices of its two components. The covariance matrix of Ujj^ is Rjj^ , and the covariance 

matrix of Udep is Rdep , therefore the covariance matrix of U can be expressed: 
Ru ~ roRoEP + (1 - ro)RiND . where ro is a scalar in the range 0<ro<L 

The optimal value for ro may then be found using a single scalar feedback loop, as described 
20 below. 

Referring again to Fig. 4, the independent noise covariance matrix and the 

dependent noise covariance matrix Rdep are established by the independent noise covariance 
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matrix sub-module 44 and the dependent noise covariance matrix sub-module 42, respectively. 
The independent noise covariance matrix is preferably calculated during the manufacture 

of the mobile device receiver and stored within a memory device accessible by the AGMF 
Weight Determination module 34. For instance, if the independent noise of concern is limited to 
Additive White Gaussian Noise, then the impulse response of the receiver during manufacture 
will yield the independent noise covariance matrix Rn^. The dependent noise covariance 
matrix Rdep is preferably calculated during operation of the mobile device by monitoring the 
channel impulse response h and the delay vector 3. Various methods for calculating the 
independent and dependant noise covariance matrices are known, and should be apparent to 
those skilled in the art. 

The independent noise covariance matrix Rjj^ and the dependant noise covariance 
matrix Rdep are provided as inputs to the total noise covariance matrix sub-module 40 to 
establish two components of R^ . The scalar parameter ro is established by the optimizer module 

46, and is provided as a third input to the total noise covariance matrix sub-module 40. In a 
preferred embodiment, the optimizer module 46 determines the optimal value for r© by first 
choosing an arbitrary or estimated value for r©, and then incrementing or decrementing ro until a 
feedback signal, such as the SNR of the decision variable z, reaches its peak or optimal value. 
An exemplary method for calculating the optimal value for ro using the SNR of the decision 
variable z is described below with reference to Figs. 5-7. It should be understood, however, that 
the scalar parameter ro could alternatively be calculated by monitoring some other feedback 
parameter, such as the bit error rate of the decision variable z, that indicates when the weight 
signals w are optimal. 
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The weight determination sub^module 50 receives the total noise covariance matrix Ru 
and the channel impulse response h , and calculates the weight signal vector w according to the 
equation described above. The weight signal vector w is then coupled as an input to the rake 
receiver 10, and settles to its optimal value, Wopt , as the scalar parameter r© is incremented or 
decremented by the optimizer module 46. 

Fig. 5 is a graph 60 plotting the SNR of the decision variable output z from the rake 
receiver 10, SNRz, as a function of the scalar parameter Tq. SNRz is represented along the y-axis 
of the graph 60, and the incremental values of ro, between zero (0) and one (1), are shown along 
the X-axis. The peak SNRz 62 corresponding to the optimal weight signal vector Wopt appears at 
the apex of the plotted curve 64 and corresponds to the optimal value for ro. In this embodiment, 
the range of r© is divided into eleven (11) discrete values, {0, 0.1, 0.2, 0.3,..., 1.0}, which are 
identified on the x-axis of the graph 60 by eleven points or states q that are labeled {1, 2, 3, 
11}. The eleven (11) discrete values of ro are therefore referred to hereinafter as states one (1) 
through eleven (11). In alternative embodiments, the precision of the optimizer module 46 could 
be increased or decreased by varying the number of states. 

Cross-referencing Figs. 4 and 5, the SNR module 48 calculates SNRz for the current state 
q, which the optimizer module 46 preferably stores at a memory location on the device. The 
optimizer module 46 then increments or decrements the current state q, and compares the new 
SNRz with the stored value. In one alternative embodiment, two consecutive states of the SNRz 
may be calculated simultaneously by using a Dual Decision Statistic Pilot Rake Receiver 70 as 
described below with reference to Fig. 6. In either case, if the SNRz of the new state is greater 
than the SNRz of the current state, then the optimizer module 46 sets the new state as the current 
state. This process may be repeated until the SNRz reaches its peak value 62 in order to achieve 



11 



the optimal value for Tq. This method for determining the optimal value for ro is described below 
in more detail with reference to Fig, 7. 

Fig. 6 is a block diagram of an exemplary Dual Decision Statistic Pilot Rake Receiver 70. 
This rake receiver 70 could be used, for example, to simultaneously provide two states of the 

5 decision variable output, z(q) and z(q'). This rake receiver 70 is similar to the rake receiver 10 
described above with reference to Fig. 1, except the array of correlator outputs {y(di), y(d2), 
y(dj)} are coupled to two output stages 72, 74. Each output stage 72, 74 includes a plurality of 
weight multipliers 76, 78 and an adder 80, 82. The weight multipliers 76 in one output stage 72 

p are coupled to a first vector of weight signals w(q) corresponding to a first state q, and the 

t§ weight multipliers 78 in the second output stage 74 are coupled to a second vector of weight 
signals w(q') corresponding to a second state q'. The outputs from the weight multipliers 76, 78 
are coupled to the adders 80, 82 in their respective output stages 72, 74, thus generating the two 

III States of the decision variable output, z(q) and z(q'). 

+^ Fig. 7 is a flow diagram illustrating an exemplary method for calculating the optimal 

weight signal vector Wopt . The method begins at step 92. In step 94, a current state value q is 
initialized. The initial value for q may, for example, be set zero, to some pre-selected estimate, 
or to any other value within the range of state values. Then, in step 96, a candidate state q' is 
estimated based on the current state q. The estimation for q' may, for example, be achieved by 
alternately selecting the state below and above the current state q each time the method 90 is 

20 repeated. For instance, in one pass through the method 90, the candidate state q' may be set to 
q+1, and then in the next pass the candidate state q' may be set to q-L In another embodiment, 
the current state q could be initialized at either the highest or lowest state value, and then 
incremented or decremented each time each time the method 90 is repeated. 
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In step 98, two weight signal vectors w(q) and w(q') are calculated, as described above, 
with w(q) corresponding to the current state and w(q*) corresponding to the candidate state. In 
an embodiment utilizing the Dual Decision Statistic Pilot Rake Receiver 70 or a similar rake 
receiver, the two weight vectors w(q) and w(q') may be calculated simultaneously. In other 
embodiments, however, the weight vectors w(q) and w(q*) may be calculated in succession 
using a rake receiver with a single output stage, such as the rake receiver described above with 
reference to Fig. 1. 

Once the weight vectors w(q) and w(q') have been calculated and applied to a rake 
receiver, the current and candidate decision statistics z(q) and z(q*) are sampled from the output 
of the rake receiver. Similar to the weight vector calculation described in step 98, the decision 
statistics z(q) and z(q') may be sampled simultaneously or successively, depending upon the type 
of rake receiver. Once a sufficient number of samples have been calculated such that the SNRz 
may be calculated with statistical significance (step 102), the current and candidate SNRz(q) and 
SNRz(q') are calculated in step 104. 

If the current SNRz(q) is greater than the candidate SNRz(q') (step 106), then the current 
state q is set to the candidate state q' (step 108). Then, in step 110, the current state q is stored, 
and the method repeats at step 96. 

This written description uses examples to disclose the invention, including the best mode, 
and also to enable any person skilled in the art to make and use the invention. The patentable 
scope of the invention is defined by the claims, and may include other examples that occur to 
those skilled in the art. 
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